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ABSTRACT 

Thermal emission from the accretion disc around a black hole can be polarized, due to 
Thomson scattering in a disc atmosphere. In Newtonian space, the polarization angle 
must be either parallel or perpendicular to the projection of the disc axis on the sky. 
As first pointed out by Stark and Connors in 1977, General Relativity effects strongly 
modify the polarization properties of the thermal radiation as observed at infinity. 
Among these effects, the rotation of the polarization angle with energy is particularly 
useful as a diagnostic tool. 

In this paper, we extend the Stark and Connors calculations by including the spec- 
tral hardening factor, several values of the optical depth of the scattering atmosphere 
and rendering the results to the expected performances of planned X-ray polarimeters. 
In particular, to assess the perspectives for the next generation of X-ray polarimeters, 
we consider the expected sensitivity of the detectors onboard the planned POLARIX 
and IXO missions. We assume the two cases of a Schwarzschild and an extreme Kerr 
black hole with a standard thin disc and a scattering atmosphere. We compute the 
expected polarization degree and the angle as functions of the energy as they could be 
measured for different inclinations of the observer, optical thickness of the atmosphere 
and different values of the black hole spin. We assume the thermal emission dominates 
the X-ray band. Using the flux level of the microquasar GRS 1915+105 in the thermal 
state, we calculate the observed polarization. 

Key words: polarization - relativity - X-rays: binaries - instrumentation: polarime- 
ters 



1 INTRODUCTION 

The imprints of General Relativity (GR) effects on the ra- 
diation emitted in the inner regions of the accretion disc 
around the black hole in Active Galactic Nuclei (AGN) and 
Galactic black hole systems (GBHS) have been the subject 
of intense observational efforts in the recent past. The iron 
line spectroscopy has been so far the most used and success- 
ful technique to study GR effects and to estimate the black 
hole spin in both classes of objects (see iMillerl 1200 7l for a 
recent review and references therein). Recently, continuum 
measurements have a lso been shown to be a promising tech- 
nique for GBHS (e.g. lLiu et aT]|2008l ). 

X-ray polarimetry is also a potentially very power- 
ful technique in this respect, as shown by iMattl (|l993l ) 
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and Dovciak, Ka r as fc Matt (l2004al) for AGN. F or GBHS, 
IStark fc Connors! (| 19771) ;IConnors fc Starkl (|l977l) and Con- 
nors, Stark & Piran (|l980l) discussed the GR effects on the 
polarization properties of the thermal emission of the accre- 
tion disc. They showed that a strong variation of the po- 
larization angle and degree with energy, due to the radial 
dependence of the temperature, is induced. In this paper we 
reconsider and extend the Stark and Connors calculations, 
mainly we consider the scattering atmosphere with different 
optical depths. 

We assume the Kerr metric for the gravitational field, a 
standard thin disc with a scattering atmosphere for the ac- 
cretion flow and Thomson scattering in the disc atmosphere 
as the origin of polarization. The resulting polarization an- 
gle and degree will be shown as functions of the energy for 
different inclinations of the observer, optical thickness of the 
atmosphere and black hole spins. 
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Figure 1. The dependence of the local polarization degree on the 
cosine of the emission angle. Curves for different optical depths 
r are shown. The results for r > 5 do not differ much from the 
case t = oo. 
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Figure 2. The multicolour black body spectra (energy flux f a E) 
for the extreme Kerr (solid line) and Schwarzschild (dashed line) 
black hole for the three observer inclinations of 30° , 60° and 85° 
and for the semi-infinite atmosphere. 



The advent of a n ew generation of X-ray polarime- 
ters (|Costa et alj l200lh makes the observation of the ef- 
fects discussed in this paper feasible. We will discuss the 
observational perspectives for a bright GBHS, the well- 
known high inclination microquasar GRS 1915+105, for 
both a relatively sm all polarimetric mission like POLARIX 



Costa et al. 20061) and for the polarimeter aboard IXO 
Bellazzini et al.ll2006bl ). 



2 ASSUMPTIONS OF THE MODEL 

We assume a Keplerian, geometrically thin and o pticall y 
thick disc around a black hole (Shakura & Svunvaev ll973l ). 
A full Gene ral Relativity treatment is adopted (Page & 
Thorne[l974|). We will show results for both a Schwarzschild 
(a = 0) and an extreme Kerr (a = 1) black hole, where a is 
the dimensionless angular momentum per unit mass of the 
black hole. At each radius the disc is assumed to emit as 
a black body, modified by Thomson scattering (see below). 
The radial dependence of the temperature is given by the 
Novikov-Thorne relativistic expression to which we add the 
hardening factor f co \ correction, 

1/4 



T(r, a) w 741 / col 
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M Q /y 
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where R(r,a) = r~ s/4 [£(r,a)/C(r,a)] 1/4 and C(r,a) and 
C(r,a) are function s of t he radius and black hole spin (see 
Novikov & Thorne Il973l for the details). We suppose that 
there is zero torque at the inner edge of the disc, which is 
assumed to be coincident with the innermost stable circular 
orbit, i.e., nsco = 6r g (for Schwarzschild black hole) and 
nsco — lr s (for an extreme Kerr black hole) with r g = 
GM/c 2 . 

We also assume that photons are scattered in the atmo- 
sphere of the disc and thus the observed radiation is polar- 
ized. We assume Thomson scattering and no absorption in 
the atmosphere. The resulting polarization in the local ref- 
erence frame of the disc atmos phere is computed using the 
Monte Carlo code STOKES (|Goosmann fc Gaskelll [20071 ) 
for finite optical depths; for infinite optical depths we adop t 
Chandrasekhar's analytical formula (Chandrasekhar 1 1960h . 



In Fig. [T] the polarization degree we obtain for the local disc 
emission is plotted versus the cosine of the emission angle, 
[i e , and for different vertical optical Thomson depths, r, of 
the scattering atmosphere. For r < 1 the polarization vector 
of the emission is aligned with the projected symmetry axis 
of the accretion disc. At higher optical depths, the direc- 
tion of the polarization vector switches to a perpendicular 
one with respect to the symmetry axis. The intermediate 
case with r = 2 has parallel polarization for /j, e > 0.3 and 
perpendicular polarization otherwise. 

The effect of hardening of the energy of photons due to 
scattering is modelled by a hardening factor, whic h incre ases 
the effective temperature (Shimura & Takahara 1 19951) . In 
Fig. [5] we show the multicolour black-body energy flux in- 
tegrated over the whole disc as it would be measured by a 
distant observer. The normalization of the flux scales with 
the cosine of the inclination angle. For the extreme Kerr 
black hole, the disc can extend closer to the black hole hori- 
zon and therefore reach higher temperatures. Therefore, the 
emission is stronger and harder in this case. Here, and in all 
computations we assume the mass of the central black hole 
M. = 14M Q , the accretion rate M = 1.4 X 10 18 g/s and the 
hardening factor / = 1.7. With this choice, the maximum 
temperature in the disc is 0.389 keV at the radius 10 r g in 
the Schwarzschild case and 1.985 keV at lr g in the extreme 
Kerr case. 

The assumption of a pure scattering atmosphere is the 
most critical one as far as the polarization is concerned. If 
the absorption opacity of the atmosphere is not negligible, 
the polarization degree is reduced and the polarization an- 
gle, for large optical depths, can even switch from perpen- 
dicular (as for a pure scattering atmosphere) to parallel to 
the p roject ion of the axis on the sky (e.g. Laor, Netzer & 
Piran 1 19901 ; Matt, Fabian & Ross 1 19931 ). Our assumption is 
therefore optimistic; a more realistic calculation would need 
a detailed modelling of the disc atmosphere. We stress here 
that, whatever the absorption opacity is, in the Newtonian 
case (i.e. when the effects of special and general relativity 
are neglected) the polarization, for symmetry reasons, must 
be either parallel or perpendicular to the projection of the 
disc axis onto the sky. The continuous variation of the po- 
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larization angle with energy (see the next section) is a clear 
signature of the GR effects. 

Once the photons leave the atmosphere the polarization 
vector can be rotated due to the strong gravity of the central 
black hole. Th e emis sion is amplified by the tran sfer fun ction 
(Cunningham 1 19751 ; Dovciak, Karas & Yaqoob l2004bh and 
the energy of photons is shifted by the gravitational and 
Doppler effects. 

We assume that scattering occurs in a non-expanding 
atmosphere, so we use the transfer function computed for 
the equatorial plane and we do not include the additional 
rotation of the polarization angle that co uld be caused by 
the bulk outflow velocity (see Beloborodov ll998T l . 

Also one can expect that magne tic fields of the disc at- 
mosphere (Silantev & Gnedin 120081 ) will interfere the ob- 
served polarization. However, the energy-resolved depen- 
dency of the thermal component exhibits a rather unique 
profile which distinguishes the gravitational effects from 
other influences. 

As the most promising polarimeters presently designed 
work in X-rays above a few keV, in this paper we limit our- 
selves to investigate the case of a stellar mass black hole (for 
AGN, the thermal emission is dominant in the UV band). 



3 NUMERICAL IMPLEMENTATION 

The main current application of our code is within the xspec 
data-fitting package, however, its applicability goes beyond 
that: the code can be compiled in a stand-alone mode, ca- 
pable of examining polarimet rical q uantities in the strong- 
gravit y regime (see Dovciak |2004| ; Dovciak et al. l2004bl ; 
l2004cr >. 

Properties of radiation are described in terms of pho- 
ton numbers. For a point-like source, the measured flux in 
the solid angle df2 is associated with the detector area 
dS a =D 2 dfl . This relation defines the distance D between 
the observer and the source. The observed photon flux per 
unit solid angle in the energy bin {E,E + AE), Af = 
dN/dtd0. o , is 



A/ = / drr 



r<t>+£±4> r^-t 
/ ^ / 

J<t> JE/g 



(E+AE)/g 



dEi fi G, 



(2) 



where G = g 2 l/-b e is the transfer function with g, I and fj, e 
being the energy shift, lensing and cosine of the emission 
angle, respectively, and 



fi = 



diV, 



dr dSj dfii dEi 



(3) 



is a local photon flux emitted from the surface of the disc. 
For a derivation of the equation ((2)1 and a more detailed de- 
scription of numer ical computations of the tra nsfer function 
see Dovciak (|2004 ) and Dovciak et al. (|2004d\ 
Let us define the specific Stokes parameters, 
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where /„, Q v , U v and V„ are Stokes parameters for light 
with frequency v (further on, we drop the index v but we 
will always consider these quantities for light of a given fre- 
quency) . 

The integrated specific Stokes parameters (per energy 



bin (E, E + AE)), i.e. Ai , Aq , Au and Av are the quan- 
tities that the observer determines from the local specific 
Stokes parameters i\, q\, u\ and v\ on the disc in the follow- 
ing way: 
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Here, dS = rdr dip is the integration element of the coordi- 
nate area on the disc and *!/ is a change of the polarization 
angle. The integration boundaries are the same as in eq. @. 
As can be seen from the definition, the first specific Stokes 
parameter is equal to the photon flux, therefore, eqs. ([2]) and 
(JSJ) are identical. 

In eqs. ©-(18} we used a law of transformation of the 
Stokes parameters by the rota tion o f axes, see eqs. (1.185) 
and (1.186) in Chandrasekhar (| 19601 ). 

An alternative customary way of expressing polariza- 
tion is in terms of degree of polarization P and two polar- 
ization angles \o and £ , defined by P Q — ^ q% + u% + v%/i , 
tan2xo = u /q , and sin2£ = v a /\/q 2 + u% + v%. However, 
these relations can be used only when we are interested in 
energy dependence of polarization. If we want to compute 
the polarization in certain broad energy range we need to 
integrate the Stokes parameters and use them instead of the 
specific Stokes parameters to get the correct result. 

Various physical effects can influence polarization of 
light as it propagates towards an observer. Here we examine 
only the influence of the gravitational field represented by 
the vacuum Kerr space-time. The change of the polarization 
angle is defined as the angle by which a vector parallelly 
transported along the light geodesic rotates with respect to 
some chosen frame. We define it in this way because in vac- 
uum the polarization vector is parallelly transported along 
the light geodesic. This angle depends on the choice of the 
local frame at the disc and at infinity. At the disc we con- 
sider the local frame co-moving with it with the a:-axis in 
the direction —d/d6 in the plane defined by the normal of 
the disc and the momentum p£ of the emitted photon. It 
is perpendicular to the momentum p£. The y axis lies in the 
plane of the disc, is perpendicular to the momentum p% and 
has direction —d/dip. At infinity we consider static frame 
attached to the observer's sky with i-axis identified with 
the impact parameter j3 (defined as positive in the direction 
—d/d6) and y-axis identified with the impact parameter —a 
(a is defined as positive in the direction d/dtp). 

The change of polarization angle \& is (see Connors et 
al. ll98dl 

(9) 



tan* 
where 
X = 
Y = 



X 



-(a — asin# )fti — (3k2 . 
(a — asin# )ft2 — /3/tj . 



(10) 
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with the black hole spin a positive when the black hole ro- 
tates counter-clockwise, i.e. in the direction d/dtp. The angle 
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Figure 3. Contour graphs of the change of the polarization angle for Schwarzschild (top) and extreme Kerr (bottom) black hole. The 
observer inclination is O = 30°, 60° and 85° (left to right). The observer is located to the top of the pictures. The innermost stable 
circular orbit is shown for the Schwarzschild case. The critical point, where the photons are emitted perpendicularly to the disc, is shown 
by a cross. The black hole rotates counter-clockwise in the Kerr case. The graphs are represented in the coordinates x = r cos tp, y = r sin ip 
in the equatorial plane where r and <p are Boyer-Lindquist coordinates. 



6 is the observer's inclination, Ki and K2 are components of 
the complex Penrose- Walker constant of pa rallel transport 
along geodesic K pw (Walker & Penrose 1 197Cll ) 

m = arp<,j -r[ap e -(r +a )PoJ/ 

-r(r 2 + a 2 )p e c .r, (12) 

«a = -rplf + r [p* - apZ]f + arplf . (13) 

Here, the polarization vector / M is 

with f/ M being the four-velocity of the co-rotating accretion 
discQ Thus it is chosen to be in the direction of the a;-axis 
in the above defined local frame on the disc. 



Below the innermost stable circular orbit we consider a free fall 
with the conserved energy and angular momentum of the matter 
on the marginally stable orbit. Note however, that we use this 
region only when computing the change of the polarization angle 
on Fig. \3\ The thermal emission is neglected in this region. 



4 POLARIZATION AT DETECTOR 

The local polarization induced by scattering (see Fig. [1} 
is denned in the local frame co-moving with the disc. The 
intrinsic polarization angle is perpendicular to the disc axis 
for large optical depths and parallel for small optical depths. 
Photons are travelling in a curved space-time, which changes 
their polarization angle as seen at infinity, eq. (0, depending 
on the particular trajectory of each photon. This results in 
a depolarizing effect, when all photons will be eventually 
summed up at infinity, see eqs. ||S)— (5). 

Due to aberration the change of the polarization an- 
gle for each photon is different at inf inity even in the 
Schwarzschild space-time (|PineaultHl977l ). In the Kerr case 
the rotation of polarization vector because of the gravita- 
tional dragging is added. The dependence of the change of 
the polarization angle on the position of the emission from 
the disc is shown in Fig. [3] It can be seen that at large dis- 
tances from the black hole this change is small. Therefore 
the depolarizing effect of the integration over this part of the 
disc will be also relatively small. On the other hand, below 
the critical point, where the light is emitted perpendicularly 
to the disc (see Fig. [3} , the change of the polarization angle 
can acquire any value. Thus the depolarizing effect of this 
region is large. Therefore for the overall polarization mea- 
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Figure 4. The energy dependence of the polarization degree (top) and polarization angle (bottom) as the observer at infinity would 
measure them for different optical depths of the atmosphere t = oo, 2.0, 1.0, 0.5 (from left to right), in the case of the Schwarzschild 
(dashed lines) and extreme Kerr (solid lines) black holes and observer inclinations of 30°, 60° and 85°. A polarization angle of 0° 
represents the direction aligned with the projected symmetry axis of the disc. 
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Figure 5. The dependence of the polarization degree (top) and polarization angle (bottom) at infinity on cosine of the observer's 
inclination for different optical depths of the atmosphere r = oo, 2.0, 1.0, 0.5 (from left to right), in the case of the Schwarzschild 
(dashed lines) and extreme Kerr (solid lines) black holes. A polarization angle of 0° represents the direction aligned with the projected 
symmetry axis of the disc. 
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sured at infinity it is important if this critical point is above 
or below the marginally stable orbit. 

It is clear that the depolarizing region is smaller for 
lower spin of the black hole. Note also that the area of the 
depolarizing region is larger for lower inclinations of the ob- 
server (the critical point moves farther away from the black 
hole). 

The decrease of the temperature with the disc radius 
implies that, the closer to the black hole the radiation is 
emitted, the harder it is. As a consequence, the effects of 
GR are larger at higher energies. The energy dependence 
of the polarization degree and angle at infinity is shown in 
Fig. [4] In the low energy limit, i.e. for E < 0.1 keV, where 
the radiation is emitted far away from the black hole, the 
polarization degree and angle are basically equal to the local 
ones for the emission angle equal to the inclination of the 
observer (compare with the Fig. [1} ■ In general, the polariza- 
tion is highest for these energies because, as said above, GR 
effects tend to depolarize the radiation emitted from regions 
closer to the central black hole. 

For high energies (above 10 keV), the polarization 
degree also increases. The polarization here is influenced 
mainly by that region of the disc, where the transfer func- 
tion is the largest and the temperature is the highest. This 
area is not very large and thus the span of the change of the 
polarization angle is small. Therefore the depolarizing effect 
is not very large. Note, however, that for the same reason 
(small area) the flux for this interval of energy is also very 
small (see Fig. [2| . 

The dependence of the polarization degree and angle 
on the inclination of the observer is shown in Fig. [5] Both 
quantities are integrated over the whole energy range (0.01- 
100 keV). The polarization degree in most cases increases 
with the inclination angle of the observer. The polarization 
angle does not change much with the observer's inclination 
and optical depth, apart from the highest inclinations. Its 
value depends, quite strongly, on the spin of the black hole. 



5 OBSERVATIONAL PERSPECTIVE 

The development of X-ray polarimeters based on the pho- 
toelectric effect has renewed the interest in the polarimetry 
of X-ray astronomical sources. Until now, the only p ositiv e 
detection is that of the Crab Nebula (Weiss kopf e t al. ll978f ). 
The Gas Pixel Det ector (GPD, Costa et al. 1200 ll ; Bellazzini 
et al. |2006al , l2007t) is one of the most advanced project in 
this field. Its current design allows a large increase of perfor- 
mances with respect to the instruments built so far, with a 
sensitivity peaking at about 3 keV and a response between 
2 and 10 keV. 

The GPD has been proposed as a focal plane instru- 
ment, coupled with an X-ray optics, for several satellite mis- 
sions. Here we discuss two feasible scenarios: 

• a relatively small, low-cost pathfinder mission, possibly 
dedicated to the X-ray polarimetry, which could perform 
the study of Galactic and bright extragalactic sources. Long 
pointed observations are possible in this scenario and po- 
larimetry at the level of ~1% could be reached in a few days 
of observation, depending on the flux of the source; 

• a large observatory, which could be dedicated to the 
study of the faint extragalactic sources and to the detailed 



investigation of the most interesting ones, following the re- 
sults from the pathfinder missions. In this scenario, the lower 
time dedicated to the X-ray polarimetry will be balanced by 
the larger area of the optics. 



We use two examples from the two scenarios above to 
verify if the effects presented in Sec. [3] will be measurable 
by the next missions with the GPD on-board. In the for- 
mer profile, various missions are currently under study and 
hence there are concrete possibilities that at least one mis- 
sion with an X-ray polarimeter based on the GPD will be 
launched in a few years. One of the options is POLARIX 
(Costa et al. 120061 ) . an Italian mission dedicated to the X-ray 
polarimetry. Its design is now under Phase A study by the 
Italian Space Agency and, if eventually selected, POLARIX 
will be launched between 2012 and 2014. The baseline as- 
sumes three optics and even an improved version with six 
telescopes is considered. As a benchmark for the pathfinder 
missions, we considered this improved design of POLARIX, 
with the warning that, for the baseline version, the same 
results can be achieved by doubling the observing time. 

With regard to the large mission scenario, the GPD was 
also inse rted in the focal plane of the XEUS mission (Bellazz- 
ini et al. l2006bl ), which was proposed to ESA for the Cosmic 
Vision 2015-2025 and has passed the first selection step. It 
is now considered as a possible focal plane instrument for 
the International X-ray Observatory (IXO) mission, arising 
from the merging of the former XEUS and Constellation-X 
missions. Because the details of the IXO mission have not 
been fixed yet, we use the XEUS design (e.g. Arnaud et al. 
120081 ) as benchmark for the large observatory scenario. 

We choose the prominent Galactic black hole system 
and microquasar GRS 1915+105 to see how the next mis- 
sions could investigate the GR effects on the degree and 
angle of polarization presented in Sec. U This source is 
very bright in the thermal state and the inclination of 
70° (Mirabel & Rodriguez 1 19941) favours a high degree 
of polarization (see Fig. 3}. We consider the spect rum i n 
a low-luminosity thermal state (McClintock et al. 120061 ). 
L < 0.3L c dd, to minimize contamination by a nonthermal 
component of emission. The dependence of the polarization 
with energy is that calculated in the Sec.Uwith r = 1.0 and 
80 = 70°. The presented results are computed for the hard- 
ening factor / = 1.7, mass of t he bla ck hole M. = 14Mq 
(Greiner, Cuby & McCaughrean |200 ll ) and for th e accretion 
rate M = 1.4 X 10 18 g/s (McClintock et al. l2006l ). 

We assume the performances of the current prototype 
of the GPD simulated by means of a Monte Carlo software, 
which have been ess ential ly confirmed by recent measure- 
ments (Muleri et al. 120081 ) . The measured degree and an- 
gle of polarization are derived from convolving the response 
of the GPD with the expected degree of polarization and 
adding a Poisson background. 

The results for the pathfinder scenario are shown in 
Fig. [5] for an observation of 500 ks, while those for a large 
mission are reported in Fig. [7] for an observation lasting 
100 ks. The variation with energy is easily observed for 
both degree and angle of polarization. The Kerr and the 
Schwarzschild cases can be clearly discriminated (the errors 
shown are at the 3-ct level). 
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Figure 6. Energy dependence of the degree (left) and angle (right) of polarization, measured with an observation of 500 ks by the 
pathfinder mission. The solid and dashed curves refer to the case of a Kerr and Schwarzschild black hole respectively, with t = 1.0 and 
6*o = 70° . Errors are quoted at the 3-cr level. Note that we set the angle of polarization to zero at low energy; the actual offset of the 
observed angle depends on the inclination of the accretion disk, which is unknown at this stage. 




Figure 7. The same as the Fig. [6] but for the large mission scenario and an observation of 100 ks. 



6 SUMMARY 

We have investigated the GR effects on the polarization 
properties of the thermal radiation emitted by an accretion 
disc around a stellar mass black hole. The polarization at 
infinity is changed from its local value due to strong gravity 
effects and fast orbital motion of the disc close to the central 
black hole. The difference is particularly dramatic especially 
in the region below the critical point, where the polarization 
angle swings over the entire range of 180 degrees. 

The most interesting interval of energy is from 0.1 to 
10 keV. Below this energy band the emission comes from far 
away from the disc centre, while above this range the flux 
rapidly decreases. Note, however, that for different maxi- 
mum temperatures of the disc the value of the upper bound- 
ary of this energy interval will change. All studied energy 
dependences will be shifted to a higher energy for a hotter 
disc. The maximum temperature is a function of the black 
hole mass, accretion rate and hardening factor as can be 
seen from the eq. JTJ. 

The polarization degree is the highest for the lowest 
studied optical depth of the disc's atmosphere and in most 
cases it grows with the observer's inclination. The polar- 
ization angle does not change much with the inclination of 



the observer apart from the highest inclinations. Instead, it 
depends quite strongly on the spin of the central black hole. 

The effects of strong gravity on the polarization state 
of emergent radiation will be accessible to the next missions 
with an X-ray polarimeter based on the Gas Pixel Detector. 
We presented the case of the pathfinder mission, feasible in 
a few years, and that of a large mission. In both scenarios, 
the energy dependence of the degree and of the angle of 
polarization will clearly discriminate between the Kerr and 
the Schwarzschild black holes cases. Even in the pathfinder 
scenario, the X-ray polarimetry will be a powerful probe of 
the metric around black holes. 
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